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Abstract. - Based on dynamical density functional theory (DDFT) we consider a non-equilibrium 
system of interacting colloidal particles driven by a constant tilting force through a periodic, 
symmetric "washboard" potential. We demonstrate that, despite of pronounced spatio-temporal 
correlations, the particle current can be reversed by adding suitable feedback control terms to the 
DDFT equation of motion. We explore two distinct control protocols with time delay, focussing 
on either the particle positions or the density profile. Our study shows that the DDFT is an 
appropriate framework to implement time-delayed feedback control strategies widely used in other 
fields of nonlinear physics. 



Introduction. — Transport phenomena of Brownian 
particles in complex geometries are a topic receiving in- 
tense and continuos attention since decades [TH3]. A 
large number of studies has been devoted to transport in 
structured ID systems such as colloids or biomolecules in 
microchannels [4], colloids in optical potentials 016], or 
cold atoms in optical lattices [7]. Theoretical studies of 
such systems have predicted spectacular effects such as 
ratchet mechanisms in systems with asymmetric spatial 
potential [I], giant diffusion [SUSHIS] and dispersionless 
transport [10] in symmetric systems under constant ex- 
ternal bias ("tilted washboards"), and the negative mo- 
bility effect 01[TT]. Many of these effects have also been 
observed experimentally (see, e.g., [SIIIIIIHIIH]), often in- 
volving colloidal systems. A related topic is how these 
non-equilibrium phenomena can be manipulated by con- 
trol forces |13] . Particularly promising are feedback control 
schemes which depend on the state of the system. A spe- 
cial case is the time-delayed feedback control method sug- 
gested by Pyragas [H| , where the control term involves the 
difference between an output variable (the control target) 
at time t and its value at time t — r, with t being the de- 
lay time. This method is particularly suitable to stabilize 
certain, otherwise unstable (periodic) states. Moreover, a 
time delay naturally occurs in experiments involving feed- 
back control due to the lag between the collection of infor- 
mation and the feedback. Indeed, time-delayed feedback 
control is nowadays used in a broad variety of non-linear 
systems such as lasers, neural dynamics, and excitable 



macroscopic media [15]. In the area of ID transport, the 
method has already been applied, on a theoretical level, 
to Brownian motors (driven by an unbiased, time-periodic 
force) [16], systems rectified by delayed correlated noise 
[17] . and flashing ratchets (involving asymmetric, time- 
dependent potentials) [13l[T8]. A first experimental real- 
ization of a feedback-controlled flashing ratchet already 
exists [in]. Very recently, feedback strategies with delay 
have also been explored, on the basis of Langevin equa- 
tions [20l[2T], as a tool to manipulate the particle current 
in tilted washboard potentials. However, despite intense 
research we are still far away from a full understanding 
of the usefulness of control schemes in transport processes 
and its potential applications in biology and nanotechnol- 
ogy 

One open, yet very important point is the role of 
particle interactions, investigations of which have only 
started recently (see, e.g., [22]). In the present letter 
we show, for the first time, that control in an interact- 
ing, driven, overdamped system can be well implemented 
within the framework of dynamical density functional the- 
ory (DDFT) [33H15]. In the DDFT the basic dynamic 
variable is the time-dependent, continuous density field, 
where the microscopic interactions determining the un- 
derlying system of discrete Brownian particles enter via a 
free energy functional. Therefore, DDFT has a bridging 
position between microscopic (Langevin-equation based) 
models, on one hand, and mesoscopic, hydrodynamic mod- 
els for transport of continuous phases (see, e.g., 26 ), on 



p-1 



K. Lichtner et al. 



■ 
0.2 



-a) 1 


<\ 


F=0 - 


... i ... ■ 


"... i 


- F=0 






z / a 



F=0.5- 



F=3 



z / o 



16 



15 

3 p 


10 

s 

& 
10 





Fig. 1: Density profiles as a function of the particle position 
at e = 4 and AT = 4. The washboard potential [7 wb (z) (with 
fc = 7r/(8cr)) is indicated by the dotted lines. Top: profiles for 
non-driven systems (F = 0) at a) U — 15 and b) U = 3. In 
a) we have included a profile for e — 1 (dashed line). Bottom: 
profiles for driven systems at c) F — 0.5 and d) F = 3 and times 
t = (solid), t = tb (dashed) and t = IOtb (dot-dashed). The 
total length of the channel is L = 60<r. 



the other hand; a connection, which is also highlighted by 
a recent derivation of the DDFT via projector operator 
techniques [27]. In the last years, DDFT has been applied 
to a variety of driven systems such as colloids in unsta- 
ble traps [35] and sedimenting colloids [35] . Moreover, in 
a very recent work [30], DDFT is employed to study at- 
tracting colloidal particles in ID (time-dependent) ratchet 
potentials. Here we apply the method to a system of soft, 
repulsive colloids in a static, tilted washboard potential. 
We focus on the impact of time-delayed feedback control. 
Indeed, we find that the current in our strongly corre- 
lated system can be efficiently controlled via time-delayed 
feedback schemes focussing either on the average particle 
position or directly on the density profile. 

Model and Results. — Our model system consists 
of N colloidal particles in an one-dimensional channel of 
length L. Particles at positions z% and zi in the channel in- 
teract via a repulsive Gaussian potential V GCM (zi — zi) = 
£oexp[— {z\ — z 2 ) 2 /a 2 } (with e > 0), a typical coarse- 
grained potential modeling a wide class of soft, partially 
penetrable macroparticles (e.g., polymer coils) with effec- 
tive (gyration) radius a [31, 32 . The particles are subject 
to a spatially periodic, symmetric ("washboard") poten- 
tial J7 wb (z) = Uq sin 2 (fcz), where k defines the wavelength 
and Uq is the amplitude. 

Typical distributions of the system in the absence of 
an external drive are illustrated in fig. [T^, and fig. [Tp. 
Specifically, we consider the space-dependent one-particle 
density p(z) = (J2iLi $ ( z " z i)) (with (• • •} being an en- 
semble average) at two values of the interaction strength 
e = e^/k^T, where fcs and T are Boltzmann's constant 
and temperature, respectively. The particle number is 



tained numerically using the DDFT formalism described 
below, with the initial configuration (t = 0) being a sin- 
gle density peak centered in one potential well. The data 
in fig. [T^, pertain to a washboard potential of strength 
U = Un/k^T = 15. Here, thermal fluctuations are too 
weak to allow for barrier crossing within the observation 
time and, consequently, the particles remain localized in 
the potential well. This is true both for strong interactions 
(e = 4) and in a small-coupling case (e = 1), where the 
narrower peak of p(z) suggests an even more pronounced 
localization. Figure [Tp shows, for comparison, a profile 
corresponding to U = 3 (and e = 4, F = 0), yielding 
a fluid-like situation with the density becoming non-zero 
everywhere. Within our calculations, we did not iden- 
tify a threshold value of U separating the two regimes, 
consistent with the continuous dependence of the equilib- 
rium mobility on U in the corresponding (exactly solvable) 
single- particle problem [3J. 

From now on we focus on the more interesting localized 
situation depicted in fig. Here, an effective motion of 
the particles can be induced by a constant tilting force 
pbias _ p o ~ ( w }t n 2 being the unit vector in z-direction), 
corresponding to a linear potential J7 blas (z) = — Fqz. We 
choose Fq > such that the particles move preferentially 
to the right. The non-equilibrium dynamics of the driven 
system is investigated via dynamical density functional 
theory (DDFT) p3r[231[27] . where the central quantity is 
the time-dependent, one-particle density of the particles, 
p(z,t). By construction, DDFT assumes the dynamics to 
be overdamped, i.e., inertial effects are neglected. The 
exact Smoluchowski equation for p(z, t) is then approxi- 
mated such that non-equilibrium two-particle correlations 
at time t are set to those of an equilibrium system with 
density p(z, t). Neglecting, moreover, hydrodynamic inter- 
actions, one obtains the key DDFT equation 23 4251127] 
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fixed to N — J_ L , 2 dzp(z, t) = 4. The data have been ob- 



where T is a mobility coefficient (i.e., T = D^/k^T 
with Dq being the short-term diffusion coefficient), 
and F[p] is a free energy density functional deter- 
mining the effective "current" j = — TpV (SF[p]/Sp). 
Specifically, F = F id + F int + F cxt where F id [p] = 
k B Tjdzp(z,t) [In (A 3 p(z,t)) - l] is the ideal part 
(with A being the thermal wavelength), F cxt = 
J dzp{z,t) (U wh (z) + C/ bias (2)) is the external field con- 
tribution, and J 7111 ' accounts for the colloidal interactions. 
Here we employ the mean-field (MF) approximation, that 
is, J** = (l/2)/dz 1 /dz 2 p(z 1 ,i)l/ GCM (|z 1 - z 2 \)p(z 2 ,t). 
Due to the penetrable nature of the Gaussian potential 
(which allows an, in principle, infinite number of neigh- 
bors) the MF approximation is known to become quasi- 
exact in the high-density limit and yields reliable results 
even at low and moderate densities [32j . 

The impact of the external drive on the density profile is 
illustrated in fig.[TJ; and fig.QJl, where we consider two val- 
ues of the driving strength F = Fq<t /k^T '. The numerical 
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calculations have been performed in channels with closed 
boundaries, such that any drift motion is, strictly speak- 
ing, of transient character. However, by choosing channels 
of sufficient lengths (L > 60a) we have ensured that all 
shown profiles are free of boundary effects. 

In the small-drive situation depicted in fig. [Hi (F = 0.5), 
the peak barely moves within the time range considered, 
indicating that the probability for the particles to jump 
over the barriers is still small. This changes at F = 3 
(fig- El) where the drive causes a shift of the entire density 
peak to the right (here and in the following, time is mea- 
sured in units of the Brownian timescale tb = er 2 / (r/ce? 1 ), 
which is of the order of 10 _9 s for typical colloids). At 
the same time, the distribution broadens in the sense 
that now several potential wells are (potentially) occu- 
pied by particles. The influence of F is also reflected by 
the time-dependence of the mean-squared displacements 
(MSD), w(t) = iV _1 E^i<(^(*) " ^(0)) 2 )- Since, within 
the DDFT formalism, we do not have direct access to the 
particle positions, we obtain the MSD rather via the re- 
lation w(t) = dzz 2 G s (z,t). Here, G s (z,t) is the self- 
part of the van Hove correlation function measuring the 
probability that a particle moves over a distance z dur- 
ing time t. We calculate this function, as well as its dis- 
tinct counterpart Gd(z,t) measuring the time dependence 
of two-particle correlations, within the DDFT formalism 
via the test particle method [33j . Some numerical results 
for the resulting MSD are given in fig. [5^; correspond- 
ing diffusion constants D = lim^oo w(t)/2t are plotted in 
fig- Wp- Within the time range considered, we find diminu- 
tive values of D for driving strength F J$ 0.5, consistent 
with the behavior of the density profile in fig. Hfc. Upon 
further increase of F, the MSD at intermediate times first 
displays sub-diffusion (where w(t) increases slower than 
linearly with t) and then superdiffusion characterized by 
much faster growth. In the long-time limit, the systems 
reach true diffusive behavior with w(t) oc t as expected 
from the Brownian friction incorporated in the DDFT. Ex- 
emplary data for the functions G s (z, t) [at F = 3] are plot- 
ted in fig. 12;. The appearance of several peaks in G B (z,t) 
reflects the existence of several "populations" of particles, 
one consisting of particles that remain in their original po- 
tential valley and the other ones consisting of moving par- 
ticles. Clearly, the self- van Hove function strongly deviates 
from the gaussian behavior expected in a " normal" , liquid- 
like diffusing system. Coming back to fig. [2p we note that 
the maximum value of F considered here is still smaller 
than the "critical" force ^ crit = max(dU wh /dz) w 5.9, 
beyond which the potential barriers are eliminated, and 
a single particle can slide freely. Therefore, we do not 
see the maximum and subsequent decrease of D expected 
in the vicinity of F cr ;t [8j[9] (indeed, investigation of that 
range is hindered by the finite length of our system which 
eventually yields boundary effects). We recall in this con- 
text that the particles in our system are interacting such 
that one may expect deviations of the behavior of D from 
the uncorrelated case .22 . Indeed, as indicated in fig. [2b, 
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Fig. 2: a) Mean-square displacement as function of time at 
several driving strengths (e = 4, N = 4). b) Corresponding 
diffusion constants. Included are date for e = 1. Parts c) and 
d) show the self and distinct part of the van Hove function at 
F = 3, e — 4, and some characteristic times. 



a decrease of the coupling strength yields a shift of the 
curve D(F) towards larger driving strength. Thus, the 
repulsive interactions in our model support the particles 
in crossing the barriers. Finally, fig. [5J1 shows the distinct 
van Hove function, G^(z, t), as calculated by the test par- 
ticle method 33 . At t = 0, where Gd is proportional to 
the usual pair distribution function, the correlations are 
restricted to the first potential well in which the parti- 
cle where confined initially (note the correlation hole at 
t = 0). At later times Gd(z,t) develops additional peaks 
in the neighboring potential valleys at z > 0, indicating- 
pronounced spatio-temporal correlations in the driven sys- 
tem. 

We now aim at manipulating the dynamics by time- 
delayed feedback control, which involves the difference 
between an appropriate system variable (the control tar- 
get) at time t and its value at time t — r. The use of 
such closed-loop strategies in ID systems subject to tilted 
washboard potentials has been previously explored, e.g., in 
Refs. [2UH2T]. These studies investigated non-interacting 
Brownian particles by direct numerical solution of the cor- 
responding Langevin equations. Here we consider an in- 
teracting colloidal system and the dynamics is described 
by the DDFT equation 

We explore two distinct control strategies. Within 
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Fig. 3: a) Normalized first moment as a function of time for 
various control amplitudes K (see eq. (HI), b) First moment 
(solid line) and effective bias (dashed) at K — 3. c) Density 
profiles for various times (F — 3, r = tb, tstart = 1.5tb). 



the first protocol, which is similar in spirit to that in 
[20] . the control target is the first moment of the den- 
sity corresponding to the average particle position z — 
f_™ dzzp(z, t). We note that particle positions and thus, 
their average, are measurable variables in experiments of 
colloidal transport (see, e.g., refs. 0j6]). To implement 
the feedback control we supplement the external poten- 
tial entering in the free energy functional J- ext (see text 
below eq. fll])) by a linear control potential Ui(z, t; r) = 
—zFi(t;r) where 

F?(t, t) = -K Q (1 - tanh [z{t) - z(t - r)]) (2) 

is a spatially homogeneous control force. In the 
DDFT equation (JT|) this force yields a current j c = 
-Tp(z,t)F c (t,T). Note the restriction -2K < Ff < 0, 
where Ko is assumed to be positive. In the numerical cal- 
culations we switch on the control force at i s tart = 1-5tb 
(i.e., Ff(t < Wt) = 0). 

The impact of the control force Ff on the normalized 
first moment (z)t = z(t)/N is illustrated in fig. where 
we have chosen F = 3 and r — Tb- This delay time cor- 
responds roughly to an intrinsic time scale of the system. 
Indeed, as seen from fig. [2^, t = tb is within the time 
range where the crossover from sub-diffusive into diffusive 
behavior of the MSD of the uncontrolled system occurs. 

In the absence of control (K = K^ajk^T — 0) the first 
moment plotted in fig. [3ji just increases with t, reflecting 
the rightward motion expected at F = 3 (see fig. QJ1) . The 




Fig. 4: Long-time velocity as function of K for several values 
of r and e (F = 3, istart = 1.5tb). Lines are guides for the eye. 



slope of the function (z)t at large t may be interpreted 
as an average velocity v = lim t _ i . 00 (d(z) t /dt). Increasing 
K from zero, the velocity first decreases until the peak 
motion stops (i.e., the time-average of (z) t becomes con- 
stant) at K = 3. This value corresponds to a balance 
between control force and bias. Even larger control am- 
plitudes then result in a significant backward motion, i.e., 
(z) t and v become negative. 

We now consider in more detail the time-dependence of 
the feedback control. First, a significant influence on (z) t 
appears only at relatively large times t 3> t s tart (which de- 
pend, in turn, on the actual value of K). At earlier times, 
the density peak moves so quickly that z(t) 3> z(t — r) 
and the tanh- function in eq. ^ approaches 1, yielding 
Ff(t, t) w 0. With the "slow-down" of (z) t at somewhat 
later times (visible also at K = 0), the argument of the 
tanh decreases. Thus, the control sets in, yielding an ef- 
fective biasing force F cS (t) = F° + Ff(t;r) < F . The 
behavior of F eS (t) and the function (z) t is shown in fig.O) 
for the "balanced" case K = 3. One sees that the control 
becomes effective at t w 10tb- After that, the first mo- 
ment displays an oscillating behavior changing between 
small backward motion and forward motion. These oscil- 
lations yield, in turn, oscillations of F oS a /k&T between its 
maximum value, F = 3, and its minimum F — 2K = —3. 
As a consequence, the overall motion stops. Interestingly, 
the oscillations (which seem to persist in the long-time 
limit) have a period of about 5tb, that is, much larger 
than the delay time (r = r B ). Density profiles related to 
one "cycle" of (z) t and F cS (t) are plotted in fig.[3fc, where 
t/rs = 25 and 30 correspond roughly to the minimum and 
maximum of F eS , respectively (see fig.[3jD). It is seen that 
the entire particle distribution is shifted with the periodic 
changes of F cS . These shifts are accompanied by changes 
in the peak shapes, which become most asymmetric when 
F eS is minimal, i.e., the control is maximal (£/tb = 25). 
To complete the picture, we plot in fig. [4] the long-time ve- 
locity v (averaged over the oscillations of (z) t , if present) 
as function of the control amplitude. We have included 
data for different delay times r and different interaction 
(i.e., repulsion) strengths e. All systems considered dis- 
play a clear current reversal at K = 3, where the velocity 
v changes from positive to negative values irrespective of e 
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Fig. 5: a) Normalized first moment as a function of time and 
various control amplitudes K for the second control loop (see 
eq. p|) and e — 4. The inset in a) contains data for e = 1 at 
otherwise same parameters F = 3, r = tb, and t sta rt = 1.5tb. 
b) Density profiles at some characteristic times for K = 1.4. 
The washboard potential is indicated by the dotted lines. 



and t. These parameters, however, do have an impact on 
the magnitude of the velocities in the two regimes K < 3 
and K > 3. Specifically, reduction of e (at fixed r) yields 
a decrease of v as compared to the case e = 4. Similarly, i> 
decreases in magnitude when the delay time decreases (at 
fixed e) from r = tb towards r = 0.2tb- In other words, 
the time delay supports the current reversal in the parame- 
ter range considered. We also note that all of these results 
are robust, on a qualitative level, against slight changes of 
the control protocol, such as a reversal the argument of 
the transcendental function in eq. ([2]). 

The control scheme considered so far focuses on the av- 
erage particle position. However, given that the basic dy- 
namical variable in DDFT is the density profile p(z,t), it 
is interesting to briefly discuss a control loop based on that 
quantity. Specifically, we consider the potential 

U%(z, [p]) = -K z (1 - tanh [p(z, t - r) - p(z, t)]) . (3) 

The impact of this scheme on the function (z) t is shown in 
fig. EK, where the inset contains data for the weakly inter- 
acting case e = 1. For small control amplitudes (K < 1), 
the behavior of both systems, e = 4 and e = 1, is sim- 
ilar to what is observed with the previous control loop 
(eq. ([2])) in the sense that the long-time velocity decreases 
with increasing K. However, contrary to this previous 



loop, further increase of K then yields an abrupt reversal 
of the motion, that is, the particles "bounce backwards" 
in the sense that v — > — oo. Moreover, this abrupt change 
occurs at values of K much smaller than the biasing force 
(F = 3). The spatio-temporal behavior of the microscopic 
density profile p(z, t) related to the sudden reversal of mo- 
tion is illustrated in fig. \5jp. At t = 3.5tb, that is, before 
the reversal, the initial density peak centered at z = has 
extensions to the next potential wells at positive z, indicat- 
ing that the particles move rightwards. However, already 
at this time, there is a small bump at positions z < not 
present within the first control scheme (see fig. [SJ:). At 
time t — 5tb after the reversal this additional bump has 
increased significantly along with a left-ward shift of the 
center of mass. In a future publication we will analyze 
in more detail to which extent the sudden reversal and 
the associated complex behavior of p(z,t) is a true insta- 
bility. Nevertheless, already the results in fig. [5] indicate 
that a local control scheme could be more efficient, in the 
sense that much smaller perturbations K are required to 
yield current reversal, compared to control focussing on a 
space-averaged quantity. 

Concluding remarks. — In conclusion, we have 
demonstrated that the transport in an interacting, driven 
colloidal system can be efficiently manipulated by time- 
delayed feedback control. Our control goal in the present 
context was to stop or reverse the motion in a tilted wash- 
board potential below the critical driving strength -F cr it be- 
yond which the washboard becomes ineffective. We have 
shown that this goal can be achieved with different control 
schemes that involve the same delay time r (chosen equal 
to the intrinsic, Brownian time scale) but different con- 
trol targets (average particle position (z) t versus density 
profile p(z, t)). The fact that the schemes produce compa- 
rable results indicates a certain robustness of the feedback 
control method for transport phenomena in ID systems, 
consistent with earlier theoretical findings for ratchet sys- 
tems [13lfT7l - H9] and for non-interacting particles in tilted 
washboards 20 21 . In that sense, our study also supports 
the more general perspective that time-delayed feedback 
control can be extremely useful for the manipulation of 
non-linear systems |15) . 

An experimental realization of the present results seems 
possible with micron-sized colloidal particles in ID tilted 
washboard potentials created by optical (laser) fields (61121 
IT9] . For instance, the study [12] employs particles of size 
a = 1.5pm in potentials with periodicity A = 3.14pm « 
2.1cr, and similar dimensions occur in [15] . The barrier 
heights in [12] are U = 4.5—11.5, and the biasing forces are 
in the range < F ^ 100, suggesting that our parameters 
for the uncontrolled system (A = 8<r, U = 15, F = 3) are 
not unrealistic. With respect to control, we stress that in 
colloidal systems the trajectories of the individual particles 
can be monitored with a video camera [j"2"l[T§] yielding, in 
principle, both the average particle position and the full 
density field as possible control target. Indeed, the average 
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colloid position is also targeted in a recent experimental 
realization of a flashing ratchet with time-delayed feedback 
control [19] . Moreover, in these experiments the delay of 
r = 5ms (arising from the finite time required to locate the 
colloids), is much smaller than the time scale for diffusion 
(~ 300ms), consistent with our assumption r < tb. 

Beyond the actual physical behavior, our study also 
shows that delayed feedback control can be conveniently 
implemented within the framework of DDFT, a recently 
established method to describe colloidal dynamics based 
on the microscopic interactions, yet with close relations 
to mesoscopic models of continuous phases [27]. As a 
side-result of our analysis, we note that the DDFT ap- 
pears to correctly describe also the non-controlled case, 
which has already been intensively investigated by other 
methods [2J[3]. Since DDFT is an approximate theory, 
the present results for the impact of control remain to be 
tested against quasi-exact data from Brownian dynamics 
computer simulations. However, given the good perfor- 
mance of the DDFT in other contexts [28][29] we expect 
our findings to be at least qualitatively right. 

Clearly, there is a number of questions prompted by 
our study. First, one needs to explore more systematically 
the precise role of the time delay in the various control 
schemes, and of its interplay with the intrinsic timescales 
of the system. Second, given the complex behavior of non- 
controlled systems at driving strength F « F cr i t (related 
to " giant diffusion" ) [8 , 9 , 22 , it would be interesting to 
extend the present analysis (where F < -Fcrit) accordingly. 
Third, we need to better understand the role of the na- 
ture of the interactions in the context of control. Here 
we have focused on soft, repulsive interactions character- 
izing, e.g., polymer coils [3T], but also paramagnetic par- 
ticles oriented parallel by a magnetic field (in ID). What 
would happen in the presence of additional attractive con- 
tributions such as van-der-Waals- or depletion interactions 
(arising from solvent particles in the colloidal suspension)? 
Indeed, as recently shown for colloids in ID ratchet po- 
tentials (without control) [3D] , attractive interactions and 
the resulting clustering tendency strongly influences the 
transport behavior. Moreover, on a macroscopic scale, 
attractive interactions determine the wettability and as a 
consequence, the transport behavior of drops on (hetero- 
geneous) substrates; a topic which is immediately rele- 
vant in the context of microfluidics [26] and may be ac- 
cessible within the present framework after a proper re- 
interpretation of the density fields. Finally, the present 
methodology may also be used to investigate systems con- 
sisting of several species, a question relevant in the context 
of particle sorting effects [Hill], Work in these directions 
is under way. 
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